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We review recent calculations of the probability that a hard parton radiates an 
additional energy fraction AE due to scattering in spatially extended matter, and 
we discuss their application to the suppression of leading hadron spectra in heavy 
ion collisions at collider energies. 



1. Introduction 

Experiments at tlie Relativistic Heavy Ion Collider RHIC in Brookhaven 
show for the first time a significant quenching of high-pj^ leading had ron 
spectra. In particul ar, t ransverse momentum spectra for neutral pion^J^ 
and charged hadron LiJ^ are suppressed if compared to spectra in p+p col- 
lisions rescaled by the number of binary collisions. This suppression is 
most pronounced (up to a factor ~ 5) in central Au+Au collisions and 
smoothly approaches the binary scaling case with decreasing centrality. The 
azimuthal anisotropy V2{p±) of hadroproduction stays close to maximal up 
to the highest transverse momentum El Moreover, the disappearance of 
back-to-back high-pj^ hadron correlations El provides an additional indica- 
tion that final state medium effects play a decisive role in hadroproduction 
up to p±_ ~ 10 GeV. 

Parton energy loss has been proposed to account for the small nuclear 
modifica tion factor ISl the azimuthal anisotropy ^ and the disappearance 
of dijets 1^^. Several studies (see references given in indicate, how- 
ever, that in the kinematical regime relevant for RHIC {p± < 12 GeV), 
p_L-broadening, shadowing, formation time and possibly other effects con- 
tribute significantly to the high-pj^ nuclear modification as well. 

Here, we discuss the status of parton energy loss calculations and their 
comparison to data in the kinematical pj_-range probed at RHIC. 
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2. Medium-induced gluon radiation from a static medium 

Several groups | 11 | 1^ | 13 114| calculated recently the modification of the ele- 
mentary splitting processes q qg and g gg due to multiple scattering 
The inclusive energy distrib ution o f gluon radiation off an in-medium pro- 
duced parton takes the form 1^1^ 
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where gluons of transverse momentum k± < £^re included. The radia- 
tion of hard quarks or gluons differs by the Casimir factor Cji — Cp or C^, 
respectively. The properties of the medium enter eq. by the product of 
the time-dependent density n(^) of scattering centers times the strength of 
a single elastic scattering o-(r) = 2/ |a(q)P (l — e*'^ ''), where |a(q)p 
denotes the elastic high-energy cross section of a single scatterer. 

For explicit calculations, one has to approximate the path integral in 
This is done either by a saddle point approximation obtained in the 
medium-induced soft multiple scattering approximation, 
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Here, the only medium-dependent quantity is the transport coefficient 
q{S^) which characterizes the transverse momentum squared ^j? transferred 
to the projectile per mean free path A. Alternatively, one can proceed in 
the opacity expansion which amounts to expanding the path integral in 
in powers of the elastic scattering center 
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Figd compares numerical results for the gluon energy distribution ob- 
tained in these two approximations. Qualitatively, this figure can be under- 
stood by estimating to what degree initial state gluons decohere from the 
partonic projectile. Decoherence depends on the relative phase (p accumu- 
lated by the gluon due to scattering. This phase grows with the transverse 
momentum accumulated by the emitted gluon. 
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Figure 1. The gluon energy distribution in the multiple soft (LHS) and single hard 
(RHS) scattering approximation calcula ted for different values of the kinematical con- 
straints R, R respectively. Fig. from RefSSl 



For multiple soft scattering, we have 



El 




^=^AzU^L=-, (4) 



which defines the "characteristic gluon frequency" lOc — \qL^. Given the 
number iVcoh of scattering centers which add coherently in the gluon phase 
one find k\ ~ iVcoh/^^- With the coherence time of the emitted gluon, 
icoh — p- — and A^coh = — y^^^rj, one estimates for the gluon 
energy spectrum per unit pathlength 

^j(m«it) I ^jlscatt 
W ; ~ — LO — ; ~ ~ Us \ — ■ (5) 

duj dz A'coh duj dz icoh V w 

This l/^/oJ-energy dependence agrees with the small-w behaviour in FigQ 
It is cut off for uj > LOc where the phase Q is smaller than unity and the 
reduced decoherence suppresses gluon emission. 

For single hard scattering [N = 1 opacity expansion) with momentum 
transfer u, decoherence occurs if the typical gluon formation time fcoh = ^ 
is smaller than the typical distance L between the production point of the 
parton and the position of the scatterer. The relevant phase is 
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which defines the characteristic gluon energy Wc = L. The gluon energy 
spectrum per unit pathlength can be estimated in terms of the coherence 
time fcoh, 

— ~ J — - as — . (7) 
The full calculation in Fig^ agrees with this l/cj-dependence in the range 

In QCD, coUinear gluons are hard and softer gluons tend to be emit- 
ted under larger angles. As a consequence, the limitations on transverse 
momentum phase space translate into a depletion of the infrared region 
of the gluon energy distribution, seen in FigQ] The estimates (O and (O 
determine the dependence on the characteristic gluon energies lUc (wc) only. 
They do not include the constraint on transverse phase space. The latter 
is determined by the parameters'^ 

Clearly, keeping lUc (lOc) fixed and taking R (i?) to infinity amounts to the 
limit of infinite pathlength. In this limit, the projectile can accumulate 
an arbitrarily large medium-induced transverse momentum and the limit 
on transverse momentum phase space is removed. For a medium of finite 
size, however, gluon emission at angles ~ (fc±jmcd suppressed since 
the emitted gluons are sensitive to the kinematical constraint kj_ < 0{uj). 
In the multiple soft scattering approximation, this translates into 

a;2 \ujj R LU, \RJ ^ ^ 

,T(mult) 

The position of the maximum of uj — as a function of R is consistent 
with this dependence on uj, see Fig. ^ 

For the single hard scattering approximation, the corresponding esti- 
mate for the infrared cut-off lo due to transverse momentum phase space 
constraints reads 

e^.4.(^)^^^i =^ iioc^. (10) 



The position of the maximum of ^^^^j^ — in Fig.^changes oc in accor- 
dance with this estimate. We thus have a semi-quantitative understanding 
of how phase space constraints deplete the non-perturbative soft region of 
the medium-induced gluon energy distribution. This suppression of the 
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non-perturbative small-o-' contributions helps to make the calculation of 
medium-induced energy loss perturbatively stable. 

In Ref."'^'"', we compare the single hard and multiple soft scattering ap- 
proximations of in detail. In general, one finds that the gluon energy 
distribution is significantly harder in the single hard scattering approxima- 
tion. For example, the average energy loss /S.E = J dujuj ^ receives in the 
multiple soft scattering approximation a dominant contribution from the 
region lu < tOc- In contrast, the dominant contribution comes from the hard 
region a; > in the single hard scattering approximation. 

Remarkably, although both approximations emphasize different kine- 
matical regions, the observable results are quantitatively comparable if 
comparable sets of model parameters are used. To relate the model param- 
eters in both approximations, one observes that up to logarithmic accuracy 
fjF'riQL ~ qL, where HqL determines the average number of scattering 
centers within the in-medium pathlength L. This implies 

R~{ni)L)R, LOc~ {nQL)ujc- (H) 

In the multiple soft scattering approximation, one uses the transport co- 
efficient q to characterize the average transverse momentum squared per 
unit pathlength. In the opacity expansion, one specifies not only this infor- 
mation, but additionally the average number of scattering centers in which 
this momentum transfer takes place. As a consequence, the opacity expan- 
sion has one additional model parameter n^L. Numerically, we find that 
for the choice no£ = 3 in the energy density distribution, quench- 

ing weight and even the a ngul ar dependence of the average energy loss are 
quantitatively comparable!^. We found deviations from this relation for 
those quantities for which due to kinematical constraints only the region 
u) < ujc contributes in both approximations, see below. This is in particular 
the case for the comparison to RHIC data in Fig^ where oj < E < Uc- 



3. Quenching Weights 



If gluons are emitted independently by a hard parton, then the probability 
P{AE) that this parton loses an additional energy fraction AE can be 
calcualted from the normalized sum of the emission probabilities for an 
arbitrary number of n gluons which carry away a total energy AE: 



P{AE) = ^ 1 

n=0 



n „ 

n J doj. 
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The summation over arbitrarily many gluon em issions in (|12|l can be per- 
formed analytically by Laplace transformation jn general, one 
finds a discrete and a continuous part, 

P{AE)^Po6{AE)+p{AE). (13) 

The discrete weight Pq may be viewed as the probability that no additional 
gluon is emitted due to in- medium scattering and hence no medium-induced 
energy loss occurs. For finite in-medium pathlength, there is always a finite 
probability po that the projectile is not affected by the medium. For infinite 
in-medium pathlength, however, one finds lim^j^ooPo = 0. 
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Figure 2. The continuous part of the quenching weight 1131 . calculated in the multiple 
soft sca ttering limit for a hard quark (upper row) or hard gluon (lower row). Fig. from 
RefM 



In Fig. 121 we show the continuous part of the quenching weight in 
the multiple soft scattering approximation. The parton energy loss AE is 
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generally seen to increase with increasing characteristic gluon energy 10^ or 
increasing factor R. This implies that AE' grows with increasing momentum 
transfer q from the medium, and with increasing in-medium pathlength X, 
as naively expected. 

To understand the negative contributions of the quenching weight for 
small values R, one recalls that the gluon energy distribution uj^ calculated 
in is only the medium-induced modification of a radiation pattern '^'^^^ ' 
which occurs in the absence of a medium, 

UJ — UJ ; VUJ—. (14) 

aio auj dio 

By writing the corresponding probabilities (|11() in Mellin space, one shows 
that 



/•oo 

p(.tot) ^^^^ ^ p^^^ „ ^) pivac) 

Jo 



(15) 



The probability p(*°*)(Ai?) is normalized to unity and it is positive definite. 
In contrast, the medium-induced modification of this probability, P{AE), 
is a generalized probability. It can take negative values for some range in 
AE, as long as its normalization is unity, 

/>oo poo 

/ dEP{E)^po+ dEp{E) = l. (16) 
Jo Jo 

The discrete weight po is found to coincide with this normalization condi- 
tion. A qualitatively similar behaviour of the quenching weights is found 
in the single hard scattering approximation. 



4. Angular dependence of radiation probability 

The quenching weights discussed above allow to calculate the average en- 
ergy loss outside an opening angle 0, 

{AE){Q)^ / dujuj— — {lo^R^oJcL) 
J duj 

= J dEE [PiE,cOc,R^iOcL)~P{E,LO,,R^ = x^tOcL)] (17) 

This calculation is straightforward since a finite emission angle can be 
taken into account simply by reducing the kinematical constraint k± < 
UJ in ^ to k± < x'^T X = sinO. In Fig|3| we compare the resulting 
angular radiation patterns obtained in the single hard and multiple soft 
scattering approximation. In agreement with the discussion following eq. 
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above, wc find that for comparable sets of model parameters Wc, R 
and ujc, -R, uq L respectively, the multiple soft and single hard scattering 
approximations lead to a comparable angular dependence of (Ai?) (O) for 
> 10°. For smaller angles, deviations persist which one can understand 
quantitatively "'^'^ and which should be regarded as an intrinsic uncertainty 
of this calculation. 

R=R=2000 




angle (decrees) angle (degrees) 



Figure 3. The average energy loss 1171 radiated outside an angle as calculated in the 
multiple soft (lower th ree lines) and single hard (upper three lines) scattering approxi- 
mation. Fig. from Ref|ii4 



5. Application 

The quenching weight P{AE) determines the quenching factor^l 

which determines the reduction of transverse momentum leading hadron 
spectra due to medium-induced energy loss. Alternatively, the quench- 
ing weight can be used to determine medium-modified fragmentation func- 
tions 1201 

DiYf{x,Q')^ r\eP{e)^^D,/^{^^,Q'). (19) 
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Figure 4. The nuclear modification factor for Tr'^-production El compared to model cal- 
culations involving parton energy loss only. Curves present the quenching factor 1181 
and the suppression factor ll^UI obtained from medium-modified fragmentation functions. 
They are given in the limiting cases where all parent partons are either quarks (upper 
lines) or gluons (lower lines). Calculations in the multiple soft scattering approximation 
use R = 2000, ujc = 67.5 GeV, corresponding to q = 0.75 



and L = 6 fm. In the 



single hard scattering approximation, we use R = R, Uc = u>c- Fig. from Ref. 



From the se, t he reduction of leading transverse momentum spectra can 
be calculated'^ by convoluting with the perturbative hard matrix elements 
which are approximately proportional to , 



(med) 
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(20) 



In Fig01 we compare the data on the suppression of the 7r"-spectra 
to suppression factors calculated from 1)18(1 and ((2()|l . We find that both 
definitions of the suppression factors lead to quantitatively comparable re- 
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suits. Moreover, the single hard and multiple soft scattering approximation 
lead to quantitatively comparable results for suitable choices of the model 
parameters. Remarkably, we find that in the presence of kinematical con- 
straints depleting the infrared region of in Fig^ the p_L-dependence 
of the quenching factor flattens considerably. This feature is necessary to 
find in Fig^a shallow pj^-dependence which is consistent with the data. 

The application of the current calculations of parton energy loss to 
data below pj^ < 10 GeV entails significant theoretical uncertainties. On 
the one hand, there are plausible competing physics effects which may be 
relevant at p± < 10 GeV (see Ref.-'^'-' for further discussion). On the other 
hand, the high-energy "eikonal" approximation used in the derivation of 
becomes questionable if AE ^ E which is the case in FigQl This should 
motivate improved calculations of ^ for which the finite energy effects of 
the partonic projectile are taken into account. 
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